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We present a study of the decay B~ — > D*'^K*~ based on a sample of 86 million T(4S) — > BB 
decays collected with the BABAR detector at the PEP-II asymmetric-energy B Factory at SLAC. 
Wo measure the branching fraction D*°K*') = (8.3 ± l.l(stat) ± l.O(syst)) x 10^", and 

the fraction of longitudinal polarization in this decay to be Fl/^ = 0.86 ± 0.06(stat) ± 0.03(syst). 
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PACS numbers: 13.25Hw 14.40.Nd 



Following the discovery of CP violation in B-meson 
decays and the measurement of the angle (3 of the uni- 
tarity triangle Q, focus has turned towards the mea- 
surements of the angles a and 7. A precise determina- 
tion of 7 requires larger samples of B decays than are 
currently available, and is likely to be based on infor- 
mation from several decay modes. Decays of the type 
are expected to play a leading role in 
this program among these modes, those with a K* 
have distinct advantages in some of the proposed meth- 
ods [3- Decay modes into two vector mesons present 
unique opportunities due to interference between helicity 
amplitudes. It has been suggested that angular analysis 
of B- &*0K*- can yield information on 7 without 
external assumptions 0. More generally, such a study 
would be sensitive to T-violating asymmetries that probe 
physics beyond the Standard Model 5]. 

The previously available information on B^ — > 
]j*Oj^*- jg based on a sample of 15 events Here 
we report on an improved measurement of the branching 
fraction for B^ D*^K*^ , and on the first measure- 
ment of the polarization in this decay. 

Results are based on 86 million T{AS) BB decays, 
corresponding to an integrated luminosity of 79 fb^^, col- 
lected between 1999 and 2002 with the BABMi detector 
at the PEP-II B Factory at SLAC 7J. An additional 
9.4 fb^^ sample of off-resonance data, recorded at e+e^ 
center-of-mass (CM) energy 40 MeV below the T{AS) 
mass, is used to study "continuum" events, e+e" qq 
{q = u, d, s, or c). 

The BABAR detector is described elsewhere 0. Only 
detector components relevant for this analysis are sum- 
marized here. Trajectories of charged particles are mea- 
sured in a five-layer silicon vertex tracker (SVT) and a 
40-layer drift chamber (DCH) in a 1.5-T magnetic field. 
Charged particles are identified as pions or kaons us- 
ing information from a detector of internally reflected 
Cherenkov light, as well as measurements of energy loss 
in the SVT and the DCH. Photons are detected in a 
CsI(Tl) calorimeter. 

We reconstruct B~ D*^K*^ in the following modes: 
D*^ D°Tr° and L»°7; D° K-n+, K-tt+tt^, and 
_ftr~7r''"7r+7r~ ; K*~ — > Kstt'^', Ks — ^ tt+tt^; tt'^ 77 
(charged conjugate decay modes are implied throughout 
this Letter). The optimization of the event selection was 
based on studies of off-resonance data and simulated BB 
events. A key feature of the analysis is the use of a sample 
of 4500 B~ — > D*°7r~ events to validate several aspects 
of the simulation and the analysis procedure. 

We select Ks candidates from pairs of oppositely- 
charged tracks that form an invariant mass within 9 MeV 
(3cr) of the known ^ Ks mass. Each Ks candidate is 
combined with a negatively charged track, assumed to be 



a 7r~, to form a K*" candidate. We retain K*^ candi- 
dates with mass within 75 MeV of the known K*~ mass. 
The Ks vertex must be displaced by at least 3 mm from 
the K*~ vertex. This requirement rejects combinatorial 
background and is 96% efficient for real Ks decays. 

Photon candidates are constructed from calorimeter 
clusters with lateral profiles consistent with photon show- 
ers. Neutral-pion candidates are formed from pairs of 
photon candidates with invariant mass between 115 and 
150 MeV. The tt" mass resolution is 6.5 MeV. 

We select D'^ candidates in the three decay modes 
listed above. To reduce backgrounds, tracks from 
K~tt'^it'^ and — > K~'it'^it'^t:~ must have momenta 
above 150 MeV. The kaon candidate track must satisfy a 
set of kaon identification criteria that provides a rejection 
factor of about 30 against pions. The kaon identification 
efficiency averaged over all kinematically allowed momen- 
tum and polar angle is 90%. For each K~Tr'^Tr'^ 
candidate, we calculate the square of the decay amplitude 
(I^P) based on the kinematics of the decay products and 
the known properties of the Dalitz plot for this decay ^3 ■ 
We retain candidates if \A\'^ is greater than 5.5% of its 
maximum possible value. The efficiency of this require- 
ment is 76%. Finally, the measured invariant mass of D'^ 
candidates must be within 2.5a of the D° mass. 

We select candidates by combining D° candidates 
with a tt'^ or photon candidate. The tt" candidate must 
have momentum between 70 and 450 MeV in the CM 
frame. The photon candidate must have energy above 
100 MeV in the laboratory frame. We reject photon can- 
didates consistent with originating from tt'^ decay when 
paired with another photon of energy above 100 MeV. 
We require the mass difference Am = m{D*'^) — m{D'^) 
to be between 138.7 and 145.7 (130.0 and 156.0) MeV for 
D*° D°Tr° (D*° L»°7). The Am resolution is 1.1 
(6.4) MeV for the D°7r° (D°7) mode. 

At each stage in the reconstruction chain, the mea- 
surement of the momentum vector of each intermediate 
particle is improved by refitting the momenta of its decay 
products with kinematic constraints. These constraints 
are based on the known mass of the intermediate particle 
and on the fact that its decay products originate from a 
common point in space. 

Finally, we select B~ candidates by combining 1)*° 
and K*~ candidates. A B~ candidate is char- 
acterized by the energy-substituted mass mss = 

(js -f po ■ pb)^/Eq — p\ and energy difference Ai? = 

Eg — i y^, where E and p are energy and momentum, the 
asterisk denotes the CM frame, the subscripts and B 
refer to the initial T{AS) and B candidate, respectively, 
and s is the square of the CM energy. For signal events 
we expect mEs — Mb within the experimental resolution 
of about 3 MeV, where Mb is the known B~ mass. 
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We require \AE\ < 40 MeV for B candidates with 
a £)° ^ K-Tr+Tr°, and |A£;| < 27.5 MeV for all other 
modes. The AE resolution is approximately 19 MeV in 
the K~tt~^tt'^ mode and 10 MeV in the other modes. 

To reduce continuum backgrounds, we make use of the 
ratio of the second to zeroth order Fox- Wolfram mo- 
ments (i?2 < 0.4), and the angle 9^ between the thrust 
axes of the candidate and the remaining charged 
tracks and neutral clusters in the event (| cos0^\ < 0.85). 
We also make requirements on the polar angle 9g of the 
B~ candidate (|cos0|j| < 0.9), and the energy flow in 
the rest of the event. We construct a Fisher discrim- 
inant J- based on the energy flow in nine concentric 
cones around the direction of the B~ candidate [l^ . 
We require < 0.40 (0.28) for B^ candidates with a 
D*° -» D°7r° iD°-f). The energy flow, 6*5,, and 0*^ are all 
calculated in the CM frame. These requirements remove 
about 80% of the continuum backgrounds and are 79% 
(74%) efficient for signal in the i:»'^7r° (D°7) mode. 

In 16% of the events there is more than one B~ candi- 
date. In these events we retain the best candidate based 
on a algorithm that uses the measured values, known 
values, and resolutions of the mass and the mass dif- 
ference Am. 

We extract the yield of B~ D*°K*~ events from a 
binned maximum-likelihood fit to the tties distribution of 
B~ candidates. The signal distribution is parametrized 
as a Gaussian and the combinatorial background as a 
threshold function The parameters of the Gaus- 

sian are determined from the mES distribution of the 
B- D*°TT~ sample. The total signal yield is 121 ± 15 
events. The third column of Table lists the yields for 
the individual modes. Figure ^ shows the mss 

distribution of B~ candidates overlaid with the fit model. 




FIG. 1: Distributions of mES for B' D*°K*-: (a) all 
modes; (b) D*° D°tv° modes; (c) D*" modes. 

The yield from the itles fit includes contributions 
from "peaking backgrounds" , which are backgrounds 
with mES near Mb- These backgrounds arise from B 
decay modes closely related to the signal mode; e.g., 
50 _^ D*^K*~ . From a Monte Carlo simulation we es- 
timate that they contribute 6.8 ± 3.4 events to the signal 
yield, where the uncertainty reflects the limited knowl- 
edge of the branching fractions for these modes. 

The branching fraction B{B^ D*°K*~) is calcu- 



lated from 



B = 



N„ 



pk 



BB 



B 



K'- 



where N^es event yield from the mES fit; ^pk 

is the peaking background, N^b = (85.8 ± 0.9) x 10^ 
is the number of BB pairs in the data sample, Bk*- 
and Bks a-re the branching fractions for K*~ — > Ks7r~ 
and Ks tt+tt^, i is an index that runs over the six 
£)*o^£)0 jnodes considered in this analysis, e]^^p is the 
event selection efficiency, and K^.o (^S^o) is the (£>") 
branching fraction for the i-th mode. This calculation 
assumes S(T(4S') ^ B+B') = B{r{4S) B°BO). The 
Monte Carlo efficiency determination uses the value of 
the polarization reported in this Letter. 

The inputs to the branching fraction calculation are 
summarized in Table Combining the six D*^/D° 
modes, we measure a branching fraction 

B{B- D*°K*-) = (8.3±l.l(stat)±1.0(syst)) x 10""*. 

The reconstruction efhciencies for photons and charged 
tracks are understood at the level of 2.5% per photon and 
0.8% per track, based on studies of a variety of control 
samples. These are the dominant systematic uncertain- 
ties in the determination of B. The efficiencies of many 
of the analysis requirements are measured in the large 
B^ I?*°7r^ control sample. The uncertainties on B 
are listed in Table ITU 

We also compute the branching fraction B{B^ 
D*^K*~) using only events from the individual D*^ / 
modes (see Table Pi. The branching fractions measured 
in the D*^ D^j modes are somewhat lower than those 
measured in the Z?*° D'^tt'^ modes. However, in the 
B^ D*^TT^ sample, where the is reconstructed 
with identical techniques, we find agreement between 
data and expectations for the relative yields of events 
in all six modes. Thus, we ascribe the difference in the 
measured branching fractions between the modes listed 
in Table m to statistical fluctuations. 

The angular distributions for the decay chains — > 
D*°K*- followed by D*" -> D"tt° or D°-f are expressed 
in terms of three amplitudes Hq (longitudinal), and 
(transverse) , and three angles, Oe, Ok and x ^^■ The 
angle 6e {Ok) is the angle of the 13" (Ks) with respect 
to the B~ direction in the D*'^ (K*^) rest frame; x is the 
angle between the decay planes of the D*'^ and the K*^ 
in the B~ rest frame. The experimental acceptance is 
nearly independent of %. Integrating over x, the angular 
distributions reduce to 

— /'F „ cx A\Ho\^cos^0dcos^ Ok 

a COS Bod COS Ok i -^^ 

+ (|i?+|2 + |ff_|2)sin2 0csin2 0K, 
-j /'F „ cx 4|i/oPsin2 6'DCOs2 0K 

a cos Pi:) a cos w A' i oi ly i\ 

I (\TT 12 I \TT \2\l^ I 



t|_£_t(j| 0111 i/LI WO UK 

+ |2)(1 + cos2 Od) sn? Ok 



for L>*o ^ D^tt" and D*^ -> L>%, respectively. 
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TABLE I: Summary of the elements of the branching fraction calculation. Nm.Es the yield from the tues fit; Npk is the 
number of peaking background events; e\.fc is the event selection efficiency for the i-th mode; = S^*- • Bks ' B^r>*o ' ^do 
the product of branching fractions for the K* , Ks, D* , and D decays in the i-th mode. 





mode 




mode 










E Umc 


X e')(xio-^) 


B{B~ ^ D*"A'*")(xlO"*) 


All 


All 


121 ± 15 


6.8 ± 3.4 


1.6 ±0.2 


8.3 ± 1.1 ± 1.0 






All 








96 ± 12 


4.8 ± 2.4 


1.0 ±0.1 


10.2 ± 1.3 ± 1.3 


D*° 




All 








24 ±8 


2.0 ± 1.0 


0.6 ±0.1 


4.4 ± 1.7 ±0.8 












B' 




D*^ 


— > L) n 




^ R- 


7^^ 




26 ±5 


1.7 ±0.9 


(6.5 ±0.6)% 


(0.54 ±0.03)% 


8.0 ± 1.8 ±0.9 


D*° 


— > U n 




K- 


7^^ 




39 ±8 


1.7 ±0.9 


(2.1 ±0.3)% 


(1.85 ±0.15)% 


10.9 ± 2.4 ± 1.7 


D*° 






~* R- 


7^^ 




31 ±7 


1.4 ±0.7 


(2.9 ±0.4)% 


(1.06 ±0.07)% 


11.6 ± 2.6 ± 1.6 






D'' 


K- 






11 ±4 


0.1 ±0.1 


(5.7 ±0.5)% 


(0.33 ± 0.03)% 


6.8 ±2.7 ±1.0 






D" 


-> K- 


7^^ 




11 ±5 


1.7 ±0.9 


(1.9 ±0.2)% 


(1.14 ±0.12)% 


5.3 ±2.9 ±1.0 




^i3°7 


D" 


K~ 




- + - 
TT TT 


0±5 


0.2 ±0.1 


(2.5 ±0.3)% 


(0.65 ± 0.07)% 


-0.2 ±3.3 ±0.4 



TABLE IL Uncertainties on Z3(B" D*°K"). 



Source 


Uncertainty 


Statistical 


13.1% 


tt" and 7 efficiency 


6.0% 


Tracking efficiency 


4.5% 


rriES fitting assumptions 


3.8% 


Event selection criteria 


3.8% 


D**^ and D'^ branching fractions 


3.2% 


Peaking background estimates 


3.0% 


Kaon identification efficiency 


2.0% 


Ks efficiency 


1.9% 


Polarization uncertainty 


1.8% 


Monte Carlo statistics 


1.7% 


Nbb 


1.1% 


Total Systematics 


11.7% 



0.04(syst) from fits to the D*° -> D°n° and D*° 
samples, respectively. Figure |21 shows projections of the 
{Qd,Qk) distributions for the event sample. 

Combining the results from the two D*^ modes, we 
find Tl/T = 0.86 ± 0.06(stat) ± 0.03(syst). The system- 
atic uncertainty reflects the accuracy of the simulation 
(± 0.017), the uncertainty on the fraction c (± 0.017), 
the finite statistics of the simulation (± 0.010), the uncer- 
tainties related to the fit assumptions (± 0.010), and the 
uncertainty due to the assumption that the acceptance is 
independent of x (i 0.004). As a consistency check, we 
fit the Qu distribution in the B~ — > D*^-k~ sample. We 
find Fi/r = 1.00 ± 0.01, in agreement with the expecta- 
tion Fl /F = 1 from angular momentum conservation. 



The longitudinal polarization fraction F^ /F, given by 

F |i/oP±|i/+P + |i/_|2' 

is extracted from an unbinned maximum-likelihood fit 
to the (Odt^^k) distribution for events with niEs > 5.27 
GeV. The data distribution (D) is fit to the sum of distri- 
butions for longitudinally {L) and transversely (T) polar- 
ized signal events, and combinatorial background events 
(C): 

D{6D,eK) = a ■ L^Od, 6k) + b ■ T{6d, eK) + c- C{9d. Ok)- 

Here c is the fraction of combinatorial background deter- 
mined from the tues yield fit, and = 1 — a — c. Thus, 
a is the only free parameter in the fit. 

The distributions of L and T are obtained from sim- 
ulations of transverse and longitudinal decays, includ- 
ing detector acceptance effects. The distribution of C is 
estimated from data candidates in a sideband of toes 
(5.20 < TOES < 5.27 GeV). We exclude from the fit 
{9e),9k) regions where the efficiency changes rapidly: 
cosOk < —0.9 and, in the D^j mode, cosOd > 0.85. 

We find longitudinal polarization fractions F^/F = 
0.87 ± 0.07(stat) ± 0.03(syst) and 0.80 ± 0.14(stat) ± 




-1 -0.5 0.5 1 -1 -0.5 0.5 1 



00560 cos 9k 

FIG. 2: Distributions of (a) cos 6d and (b) cosOk for D*" 
D^tt". Distributions of (c) cos 9d and (d) cos 6k for Z)*" — > 
D'^j. The solid line represents the full fit model, the dashed 
line represents the transverse component, and the shaded re- 
gion represents the combinatorial background component. 
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In summary, we have measured the branching frac- 
tion for B- -> D*°K*~ to be i3 = (8.3 ± l.l(stat) ± 
l.O(syst)) X 10^'*. Our measurement is a factor of 2.5 
more precise than the previous result. It is in agreement 
with predictions based on the measured B~ D*^ p~ 
branching fraction [3] , and the value of the Cabibbo an- 
gle. We have also measured the longitudinal polarization 
fraction in this decay to be Tl/T — 0.86 ± 0.06(stat) ± 
0.03(syst). This last result is consistent with expecta- 
tions [lal based on factorization, Heavy Quark Effec- 
tive Theory, and the measurement of semileptonic B- 
decay form factors, assuming that the external specta- 
tor amplitude (6 cW*~; W*~ K*~) dominates in 
B- D*°K*-. 
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